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Abstract 

Strongly intermittent turbulence is observed in the shadow of a limiter in the Large Plasma 
Device (LAPD) at UCLA [W. Gekelman, H. Pfister, Z. Lucky, J. Bamber, D. Leneman, and J. 
Maggs, Rev. Sci. Inst. 62, 2875 (1991)]. The amplitude probability distribution function (PDF) 
of the turbulence is strongly skewed, with density depletion events (or "holes") dominant in the 
high density region and density enhancement events (or "blobs") dominant in the low density 
region. Two-dimensional cross-conditional averaging shows that the blobs are detached, outward- 
propagating filamentary structures with a clear dipolar potential while the holes appear to be part 
of a more extended turbulent structure. A statistical study of the blobs reveals a typical size of 
ten times the ion sound gyroradius and a typical velocity of one tenth the sound speed. 
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Turbulence which is intermittent, meaning patchy in space or bursty in time, is often ob- 
served in both neutral fluids and plasmas A signature of intermittency in turbulent 
measurements is a non-Gaussian amplitude probability distribution function (PDF), e.g. 
with tails caused by a greater frequency of large amplitude events. Intermittent turbulence 
is ubiquitous in the edge of magnetic confinement laboratory plasmas including tokamaks (^J , 
stellarators [4], and linear devices ,5]. The intermittency in these environments is generally 
attributed to the creation and propagation of filamentary (magnetic field aligned) struc- 
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tures jal2|- Large amplitude events in Langmuir probe [8|, beam emission spectroscopy |9|, 
and gas puff imaging measurements [l^ are therefore due to the passage of high density 
structures through the low density edge region. The outward propagation of these struc- 
tures results in significant particle transport in the edge of magnetic confinement devices 
like tokamaks In addition, recently a correlation has been found between increased in- 
termittency in the scrape-off-layer (SOL) and the density limit in tokamaks, leading some to 
suggest that catastrophic transport enhancements associated with intermittent turbulence 
may be responsible for this disruptive limit jl^ . 

Several mechanisms have been proposed for the generation of these structures inclu ding 
solitary drift wave vortices propagating avalanche-like events in the plasma core 15, 

fl6 ] . zonal-flow driven generation jl?! and interchange driven production 18]. Once the 
structures are produced and ejected into the low density region, their continued cross-field 
—ion has bee,, attributed to E x B velocity due to polarization by drift-ehargiug fl 
I19I |. In tokamaks, the drift-charging is attributed to interchange forces associated with 
magnetic field gradients and curvature 0. In linear devices, rotation [3] is usually invoked 
as an interchange force. 

In this Letter, a detailed study of turbulent structures associated with intermittent tur- 
bulence in a linear magnetized plasma with a steep density gradient is presented. In this 
study, a clear observation of both density enhancements (or "blobs") and density deple- 
tions (or "holes") is made. The production of blob and hole structures is consistent with 
an interchange-like process. However, the experiments occur in a plasma with straight, 
uniform magnetic field and without significant rotation due to the straight, vertical edge 
created by the limiter. The typical interchange forces are therefore not available, yet the 
resulting intermittency is very similar to that observed in tokamaks. It has been proposed 
that a frictional force due to neutral flows may serve as an interchange force in cases where 
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rotation or curvature are not present [20j, and this theory may explain the observations. 
Two-dimensional measurements of cross-field structure show that blobs are isolated struc- 
tures with a dipolar potential while the holes are instead part of a more extended structure. 
A statistical study of the properties of the blobs has been performed, showing that the size 
of the structures scales like the ion sound gyroradius (5 ~ 10p s ) and that the velocity of the 
objects (outward, into the limiter shadow) scales with the sound speed (v ~ C s /10). 

The experiments were performed in the upgraded Large Plasma Device (LAPD), which is 
part of the Basic Plasma Science Facility (BaPSF) 2l( at UCLA. LAPD is an 18m long, lm 
diameter cylindrical vacuum chamber, surrounded by 90 magnetic field coils. Pulsed plasmas 
(~ 10ms in duration) are created at a repetition rate of 1Hz using a Barium Oxide coated 
cathode source. Typical plasma parameters are n e ~ 1 x 10 12 cm -3 , T e ~ 5eV, Tj ~ leV, and 
B < 2kG. The working gases used in these experiments were Helium and Neon. The length of 
the LAPD plasma column can be changed by closing a hinged, electrically-floating aluminum 
plate located 10m from the cathode. In the experiments reported here, this plate is partially 
closed and used to limit the plasma column as shown in Figure^a). With the plate partially 
closed, there is no parallel (to the magnetic field) source of plasma for the region behind the 
plate and any observed plasma density must come from cross-field transport. Measurements 
were performed using radially- insertable Langmuir probes [as shown in Fig. ^a)]. Triple 
Langmuir probes with 2mm long, 0.76mm diameter cylindrical tungsten tips spaced 3mm 
apart were used to measure density, electron temperature and floating potential and to 
perform two-dimensional cross-conditional averaging measurements. A linear array of six 
single Langmuir probes (1.5mm diameter flush-mounted tantalum tips spaced 5mm apart, 
aligned with their surface normals along the magnetic field) was also used to study the 
propagation of structures. 

For these experiments, the floating-plate limiter was closed so that nearly half of the 
plasma column is blocked, and a half-cylindrical plasma is observed downstream from the 
limiter. Steep gradients in plasma density are observed behind the limiter, along with very 
large amplitude (5n/n ~ 1) fluctuations. Figure ^b) shows a measurement of the radial 
profile of plasma density (n e ) and normalized root-mean-square (RMS) ion saturation current 
(J sat ) fluctuation amplitude (<5/ rms // sat ) 0.5m downstream from the limiter (He discharge, 
1.5kG field) along with an average fast Fourier transform (FFT) power spectrum of the 
fluctuations. The observed density gradient is very steep with a scale length of order a 
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few centimeters (the ion sound radius, p s ~ 0.3cm and the ion gyroradius, pi ~ 0.15cm). 
In the limiter shadow, the density profile is flat, similar to the profiles observed in the 
scrape-off-layer of tokamaks jj. The density behind the limiter is a significant fraction of 
the core density, indicative of substantial cross- field particle transport. The fast Fourier 
transform (FFT) power spectrum of the fluctuations is broadband and free from coherent 
modes. Figure ^c) shows the profile of space potential (derived from floating potential, 
also shown) in the plasma. While there is a strong feature in the potential profile near the 
location of the limiter, indicating a sheared E x B flow, there is little or no flow in the 
region behind the limiter. It is important to note that the limiter has a straight vertical 
edge. Thus, while there are vertical flows localized to the limiter edge region, there is no 
rotation and associated centrifugal force. 

Figure |2fa) shows example J sat signals measured at three spatial locations (x = 
1.3, —0.7, —2.2cm, where x = 0cm is the edge of the limiter). On the core side of the 
gradient region (x = 1.3cm), the signal is dominated by downward-going events while on 
the low-density side of the gradient region (x = —2.2cm) the signal is dominated by upward- 
going events. The upward-going events in the low density region are density enhancements 
or "blobs" , similar to those observed in the edge of many magnetic confinement devices j^J . 
The downward-going events are identified as density depletions or "holes," which are not 
typically observed in other devices, but have been reported to be seen in measurements 
on the DIII-D tokamak [9] . Generation of holes has been observed in simulations |22j and 
inward propagation of holes has been suggested as a mechanism for impurity transport 
in tokamaks |23J]. Although the raw signals clearly show the intermittency in the fluc- 
tuations, a more quantitative measurement of the intermittency can be derived from the 
fluctuation amplitude PDF. The amplitude PDF has been computed as a function of the 
amplitude normalized to the root-mean-square (RMS) fluctuation level, and is shown for 
x = 1.3, —0.7, and —2.2 cm in figure [2^b). The existence of the holes and blobs produces 
tails in the PDF at negative and positive normalized amplitude, respectively. The PDF at 
the location of peak fluctuation amplitude (x = —0.7cm) is more symmetric, although is 
still distinctly non-Gaussian. The amplitude PDF as a function of position is shown as a 
contour plot in figure Efc). From this figure it is evident that the holes are observed only 
in a narrow region of space on the high density side of the gradient (where the negative tail 
is present, 2.8 > x > 0). However the blobs are observed everywhere to the right of the 
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density gradient region (x < —0.5), as shown by the presence of a tail in the PDF at positive 
amplitude. For x < —3cm, the PDF is blob-dominated and nearly independent of position, 
suggesting that the blob structures are long lived as they propagate into the shadow of the 
limiter (intermittent signals have been observed up to 20cm beyond the limiter edge). It 
should be noted that the PDF displayed in figure |2fc) exhibits trends which are very similar 
to those observed using BES in the DIII-D tokamak (see figure 3 of reference 9). The nature 
of the intermittency and turbulent structures in these experiments is therefore quite similar 
to that observed in toroidal confinement devices, even in the absence of magnetic curvature 
or rotation. 

Figure ED^a) shows the conditional average of many (N = 16059) blob events using 
J sat signals from one tip on the linear Langmuir array (solid black line, marked "CA"). The 
amplitude threshold for triggering event selection in this case was set to twice the root- 
mean-square fluctuation amplitude. The conditional average blob event is asymmetric in 
time, with a fast rise and slow decay. This time asymmetry is consistent with observations 
in the edge of many other magnetic confinement devices Q] and with simulations j^J. Along 
with the conditional average, figure a) shows the cross-conditional average blob event as 
measured by other tips in the Langmuir array (separations shown are in the — x direction). 
The right-most tip (furthest away from the gradient region in the blob case) is used for 
event triggering. The average blob appears to travel across the probe array out into the 
low density region with a speed of ~ 942m/s, which is ~ C s /10. Figure Efb) shows the 
same conditional and cross-conditional averaging for hole events. The average hole decays 
quickly in space, but appears to propagate back into the core plasma. The interpretation 
of the Langmuir array conditional averages as being caused by propagation of structures in 
and out of the core plasma can not be conclusive without knowledge of the two-dimensional 
structure of the objects. For example, the observations could also be explained by, e.g., the 
vertical (y) propagation of structures which are tilted in the xy plane. 

In order to conclusively determine the structure and direction of propagation of the 
holes and blobs, two-dimensional cross-conditional averaging was performed using two triple 
Langmuir probes separated along the magnetic field by 60 cm. The first probe, which was 
the reference or trigger probe, was left fixed in space while the second probe was moved 
to 441 positions in a 10cm by 10cm cross-field (xy) plane centered on the position of the 
reference probe. Figure Etc) and (e) show the 2D cross-conditional average of J sat and Vf , 
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respectively, for blob events (the fixed probe was located at x — —3cm). The blob is clearly 
an isolated, detached structure, with similar extent in the two cross-field directions. The Vf 
measurement clearly shows a dipole structure of the potential associated with the blob (theO 
peak potential value is ~ 1.3V while the minimum value is ~ —1.5V, here T e ~ 5eV). The 
potential structure is consistent with E x B propagation almost entirely in the — x direction 
with an average speed of ~ 985m/s. This value is comparable to the speed measured by the 
linear probe array. Because the blob velocity is dominantly in the — x direction, the decay 
of the cross-conditional average in figure Efa) can be attributed to a finite spread in the 
blob velocity. Figure Of d) and (f) show 2D cross-conditional averages for hole events (fixed 
probe at x = 1.5cm). The hole does not appear to be an isolated structure, but instead 
is associated with turbulent structure that is extended in the vertical (y) direction. The 
potential structure does show a tendency for E x B propagation back into the core plasma, 
consistent with the linear array measurements. 

Figure Ufa) shows the measured dependence of the time width of blob events versus 
magnetic field. The inset figure shows the conditional average / sat signal for a blob event for 
three field values. As the field is decreased, the time width of the event increases. The main 
figure shows the PDF of the time width [full width at half max (FWHM)] of blob events for 
the same three magnetic field values. To calculate the PDF of the blob time width, events 
are selected from the J sat signals and the time width (AtpwHivO of each individual event is 
measured. As the magnetic field is decreased, the peak of the PDF shifts to larger time 
width and also the width of the PDF increases. The width in time of the blob event can 
be translated to a blob size using the blob velocity as measured by the linear probe array. 
Figure Efb) shows the average blob size ((pb)) compared to the ion sound gyroradius (p s ) for 
blob events measured in several different magnetic field values in Helium and one condition 
using Neon as the working gas. In addition, the average size obtained from the 1.5kG 2D 
cross-conditional average [figure El^c)] is shown. The average blob size scales with the ion 
sound gyroradius, and is approximately ten times this scale. The observation of a gyroradius 
scaling is consistent with recent theory 2^ which attributes the scaling to blob stability. It 
should also be noted that the average blob size is also comparable to the density gradient 
scale length. 

In summary, both blobs and holes are observed to be associated with strongly intermittent 
turbulence in the shadow of a limiter in LAPD. The blobs are nearly cylindrical filamentary 
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structures which have a dipolar potential, resulting in E x B propagation into the low 
density limiter shadow. A statistical study of the blobs has been performed, showing that 
the average speed of the structures is ~ C s /10 and that the average blob size is ~ 10 p s . The 
holes do not appear to be detached, propagating objects, but are instead part of a more 
extended turbulent structure, a finding which might also apply to observations of holes in 
tokamaks. These observations are made in an experiment free from traditional interchange 
forces, yet are very similar to measurements in the edge of toroidal confinement devices. 
The formation of these turbulent structures and their polarization may be linked to the 
instabilities associated with sheared flow in the gradient region. The dipolar potential of 
the blobs may also be explained by drift-charging arising from a newly-proposed force due 
to neutral flows. An initial calculation using this theory and LAPD parameters results in a 
blob velocity which agrees with the experimental observation [20| . The role of sheared flows 
and neutrals will be explored in more detail in future research. 

Discussions with G. Antar, S. Krasheninnikov, R. Moyer, A. Pigarov and D. Rudakov 
are gratefully acknowledged. The experiments reported in the paper were performed using 
the UCLA Basic Plasma Science Facility, which is funded by NSF and DOE. This work was 
supported by a DOE Fusion Energy Sciences postdoctoral fellowship and by DOE grant 
DE-FG02-02ER54688. 
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FIG. 1: (Color online) (a) A schematic of the LAPD, showing the experimental geometry, (b) 
Spatial profile of density and fluctuations behind the limiter (x = at the limiter edge) , along with 
an average fluctuation power spectrum (log-log), (c) Spatial profile of space potential and floating 
potential, as measured using a triple probe. 
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FIG. 2: (Color online) (a) Example raw 7 sa t signals at three spatial locations, (b) Amplitude PDF 
at the same three locations, (c) Contour plot of the amplitude PDF versus position. 
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FIG. 3: (Color online) Cross-conditional average of I sa t on the linear probe array for (a) blob and 
(b) hole events, showing apparent propagation of blobs out of the plasma and holes back into the 
plasma. Two dimensional cross-conditional averages of blob ((c) I sa t and (e) Vf) and hole ((d) / sa t 
and (e) Vf). All 2D conditional averages are normalized to the maximum of the absolute value of 
the average, and the color bar in (c) applies to all images. 
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FIG. 4: (Color online) (a) PDF of blob event time width, for three values of the magnetic field. 
Inset are conditional average blob events for the three magnetic field values, (b) Mean blob size 
for several conditions in He and Ne, compared to the ion sound gyroradius. 
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